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ABSTRACT 

'^here  is  increasing  interest  in  map  features  such  as  points,  lines  and 
regions  both  as  a pictoral  data  base  for  resource  management  and  as  an  aid 
to  identifying  objects  in  aerial  images.  Owing  to  the  very  large  amount  of 
data  involved,  and  the  need  to  perform  operations  on  this  data  efficiently, 
the  representation  of  such  features  is  a crucial  issue.  We  describe  a 
hierarchical  representation  of  map  features  that  consistsof  binary  trees  with 
a special  datum  at  each  node.  This  datum  is  called  a strip  and  the  tree  that 
contains  such  data  is  called  a strip  tree.  Lower  levels  in  the  tree  corresponds 
to  finer  resolution  representations  of  the  map  feature.  The  strip  tree  structure 
is  a direct  consequence  of  using  the  method  for  digitizing  lines  given  by 
[Duda  & Hart,  1973;  Turner,  1974;  Douglas  & Peucker,  1973]  and  retaining  all 
intermediate  steps.  This  representation  has  several  desirable  properties. 

For  features  which  are  well-behaved,  calculations  such  as  point-membership 
and  intersection  can  be  resolved  in  O(logn)  where  n is  the  number  of  feature 
points.  The  map  features  can  be  efficiently  encoded  and  displayed  at  various 
resol utions All  these  properties  depend  on  the  hierarchical  tree  structure 
which  al lows  primitive  operations  to  be  performed  at  the  lowest  possible 
resolution  with  great  computational  savings.  The  strip  tree  representation 
also  can  allow  parts  of  the  map  feature  to  be  accessed  sequentially.  This 
feature  is  usually  desired  when  the  map  feature  is  used  in  analyzing  images. 

The  price  paid  for  the  improved  performance  is  an  increased  storage  cost. 
This  is  approximately  4n,  where  n is  the  storage  needed  to  represent  the  xy 
coordinates. 


The  research  described  in  this  report  was  supported  partially  by  DARPA 
Grant  #N00014-78-C-0164  and  partially  by  NIH  Grant  irR23-HL-?1253-01 . 
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1.  Introduction 

We  present  a general  representation  for  polylines  (connected 
line  segments)  and  areas  (closed  polylines) . Although  this 
representation  may  have  wide  applications,  its  principal 
motivation  arose  from  the  problem  of  representing  geographical 
data  bases  of  map  features. 

A map  has  several  interesting  liinds  of  features  such  as 
contour  lines,  la)ces,  rivers,  roads,  etc.  These  can  be  roughlv 
divided  into  four  feature  classes  for  representation  in  the 
computer  T Sloan,  1978  1; 


feature 

examples  in  map  domain 

points 

towns  (large  scale  maps) 
bridges  (small  scale  maps) 

1 ines 

roads,  coastlines 

strips 

wide  roads,  rivers 

regions 

lakes,  counties 

Our  main  interest  is  in  representing  lines  and  regions.  A point 
is  such  a simple  datum  that  it  can  be  easily  treated  as  a 
primitive  in  any  representation.  Collections  of  points  from  a 
single  class  can  be  efficiently  represented  as  h-d  trees  TBently, 
1975;  Barrow  et.al.,  1977  I and  so  points  are  not  the  focus  of 
our  interest,  although  they  do  interact  with  our  representation. 
A strip  feature  is  essentially  a line  where  a locally  varying 
thickness  is  important,  examples  of  which  are  rivers  and  roads. 
As  we  shall  see,  our  representation  for  lines  will  also 
this  type  of  feature. 


encompass 
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We  regard  collections  of  these  map  features  as  a data  tase 
that  might  le  used  to  perform  the  following  tasks: 

.7ind  where  a road  intersects  a river 

•Display  a subset  of  map  features  that  appear  in  a given  map 

sector 

•Find  out  if  a given  point  is  in  a region 

-Search  an  aerial  image  near  the  edge  of  a dock  for  ships. 

A very  important  aspect  of  all  these  tasks  is  that  we  may  he 
satisfied  if  they  are  performed  at  resolution  lower  than  the 
ultimate  resolution  represented. 

Our  representation  for  lines  and  regions  consists  of  a 
binary  tree  structure  where,  in  general,  lower  levels  in  the  tree 
correspond  to  finer  resolutions.  The  tree  structure  is  a direct 
consequence  of  using  the  method  for  digitizing  lines  given  by 
fDuda  and  Hart,  1973;  Turner/  1974  ] and  retaining  all 
intermediate  steps  in  the  digitization  process.  As  an  example  of 
the  representation,  Figure  1 shews  some  roads  represented  at 
various  levels  (resolutions)  in  the  tree  structure. 

The  idea  of  representing  a line  by  sets  of  strips  was 
recognized  ty  fPeucker,  1976  ].  In  particular  he  was  able  to  find 
line  intersection  and  point  in  polygon  algorithms.  However,  the 
tree  structure  is  a vast  improvement  over  the  se'*-  organizat  icn: 
the  algorithms  are  more  efficient,  line-area  in  te  re  section  and 
area-area  intersection  and  union  can  now  he  dealth  with,  and  the 
tree  structures  are  closed  under  these  operations. 


iqure  1.  Map  features  displayed  at  various  resolutions 


using 


the  hierarchical 


structure . 
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2.  The  Strip  Tree 
2.1  Notation 

We  define  a strip  segment  L (delta)  as  the  vector  L and  the 
scalar  delta  as  shown  by  Figure  2.  The  vector  L starts  at 
(X3eq,Y3eq)  and  ends  at  (XEnd,YEnd).  We  use  S to  denote  the  set 
of  points  inscribed  by  the  rectangle  defined  by  L (delta).  Also 
we  denote  the  boundaries  of  the  rectangle  by  the  line  segments 
i*,  1-,  e + , e-  as  shown. 


Figure  2.  Definition  of  a Strip  Segment. 

A polyline  is  an  ordered  list  of  discrete  points  yO,.,.,yn 
subsets  of  which  may  be  colinear.  For  the  moment  we  require 
these  points  to  be  considered  as  connected;  later  we  will  relax 
this  condition.  We  say  a polyline  is  represented  at  resolut ion 
delta*  if  there  exists  an  ordered  sequence  of  m strip  segments 

I.^(delta)  , k=0/  ...,  m-1 

such  that 


delta  ^<delta*  )c  = C,  «..»  m 

yi  Q LJ  i i = l»  . . . , n 

If  within  a strip  segment  there  is  a point  y that  is  a member  of 
e*,  another  that  is  a member  of  e-,  and  ^-here  is  a point  v that 
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is  a m€at€r  of  1+  and  another  that  is  a member  of  1-,  then  the 
strip  segment  is  said  to  be  compact . The  compactness  property  is 
verv  important  for  some  of  the  algorithms  which  follow.  Figure  1 
shows  some  examples  for  different  deltas. 


2.2.  Digitization 


Suppose  we  have  a polyline  ?.  such  as  shewn  by  Figure  3a. 
For  any  resolution  delta  we  can  approximate  this  line  with  strip 
segments  as  follows  f Duda  S Hart,  1973;  Turner,  19741: 

Consider  the  polyline  E defined  by  (y^,yT\3  For 

each  point  y ? find  the  perpendicular  distance  d(y) 

from  y to  ?.  Denote  the  subset  of  y ? such  that  v.I>C' 

as  ?+ . P-=  E-pt.  Now  find  d*  = max  d(y)  and  d-  = max 

yeP" 

d(y).  If  (d+)  * (d-)  <delta*  then  the  polyline  is 

compactly  represented  at  resoluticn  ^ by  the  strip  tree 
consisting  of  a single  root  strip  1 ( (d  + ) e (d-) ) . If  not 
then  the  desired  strip  tree  is  obtained  tv  recursively 
applying  the  algorithm  tc  the  Ps  ^ yO ,.  . . , y+3 

,vn3  and  making  the  results  the  left  son  and 
right  son  respectively  of  the  strip  tree.  In  the  case 
of  ties  for  the  maximum  distance  d,  we  will  arbitrarily 
pick  the  point  nearest  the  mid  point  (in  arc  length)  . 


For  the  purposes  of  the  union  and  intersection  algorithms  tc 
follow  it  is  helpful  tc  think  of  the  strip  trees  as  completely 
expanded  down  to  individual  points,  even  though  these  points  may 
be  colinear.  Figure  3 shows  an  example  of  two  levels  of 


Figure  3 shows 


?a  qe  P 


recursion  cf  this  algorithm. 


Figure  3.  Steps  in  the  Digitization  Process. 

To  see  formally  that  the  convergence  is  guaranteed,  note 
that  a ? cf  k points  can  always  ke  approximated  by  a single  strip 
segment  L |k)  with  length  k assuming  eight-connectedness.  Ihus 
for  any  delta  there  must  be  a strip  tree  with  leaves  consisting 
of  no  more  than  n/delta  strip  segments  which  approximate  ?. 
Since  the  digitization  algorithm  splits  each  ? into  two  parts 
such  that  each  part  has  finite  length,  the  process  oust 
ultimately  consider  sets  of  F of  delta  points  or  less. 

2.2  Strip  Iree  definitions 

The-  binary  tree  resulting  from  the  digitization  process  is 
called  a strip  tree,  where  the  datum  at  each  node  is  a strip,  L. 
The  nodes  of  the  tree  are  initially  ordered  on  arc  length. 
(Later  we  will  see  that  when  intersection  occurs  in  two  areas 
which  are  represented  in  strip  trees,  this  property  is  sometimes 
nc^  preserved) . 

In  the  ensuing  algorithms  we  will  use  the  following 
def init ions : 

T = Evmbcl  for  a Strip  Tree  obtained  by  the  digitization 


process . 
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•5  (T ) = *:  he  points  associated  with  the  strip  at  ’■he  root  node  of 
T;  i.e.  fxl  xC  S(T)  } 

Area  (1)  = the  area  associated  with  the  strip  at  the  root  node  of 

he  measure  area  in  pixels  so  that  a strip  L(0)  still  has 
finite  area.  The  most  primiti'.'e  strip,  a single  pointy  has 
unit  acea. 

LSon  (T)  = the  left  sen  of  the  node  T 

FSon  (T)  = t he  right  son  of  the  node  T 

A node  of  the  strip  tree  is  completely  defined  by  the  seven-tuple 
(LSen,  F.Son,  Area,  XBeg  , XEnd,  XEeg,  YInd),  The  measure  Ar6a(T) 
is  tetter  fer  some  of  the  algorithms  to  follow.  Area  and  delta 
are  related  by  de Ita  s Area/ 1 | L | | . 

2.U.  vhy  Binary  Trees? 

The  polylines  can  also  be  represented  as  a tree  with  nodes 
of  more  than  two  siblings.  In  fact, nodes  could  have  different 
r.uaters  of  siblings  which  would  still  be  ordered.  Figure  4 shows 
an  example  of  the  alternate  encoding  scheme.  In  certain  cases 
this  may  be  a more  concise  representation  for  the  polyline  and 
for  all  the  algorithms  t ht  follow  we  can  extend  the  operations 
from  two  sons  to  multiple  sons.  However,  this  change  does  not 
alter  the  complexity  of  the  operations  that  we  would  like  to 
perform  and  can  be  more  inefficient  than  the  I inary  tree 


rep  re  sen  ta  *•  ion 
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Figure  4:  A portion  of  an  encoding  using  m-ary  trees. 

3.  Operations  on  Polylines 

Coaput at ional  conplexity  of  the  various  operations  is 
difficult  to  characterize,  as  it  depends  on  the  particular 
geometry  of  polylines.  If  the  polylines  are  " weli-Leha ved " , that 
is  they  are  relatively  smooth  and  do  not  self -intersect  for  acre 
than  a few  points,  then  the  algorithms  are  very  efficient.  Vhat 
'his  means  for  a particular  operation  in  terms  of  the  strip  tree 
IS  that  if  the  number  of  strips  that  must  be  examined  at  any 

level  is  constant,  then  the  complexity  of  the  operation  is  0(loq 

n)  . 

3.1.  lestinq  the  Proximity  of  a Point 

If  we  would  like  to  find  out  if  a point  is  near  a polyline, 

this  may  te  discovered  early  using  the  strip  tree.  We  can  irake 

t.nis  more  precise  Ly  exploiting  the  following  property; 

Proper t y Pi : 

A.  If  a point  z is  inside  a compact  strip  l(delta1 
then  it  can  be  at  most  2^  units  away  from  ♦•he  P. 

E.  If  a point  z is  outside  a compact  strip  l(delta) 
then  the  distance  of  the  point  from  the  ? is  bounded  by 
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3 ^ ^ 1 l.{delta)  ) + dfeita 

It  is  interesting  to  study  these  bounds  as  the  depth  in  the 
resolution  tree  increases.  Although  the  convergence  is  not 
aonotonic,  the  bounds  do  converge  to  the  actual  set- t heoretic 
distance  d^(z,P).  Now  suppose  we  want  to  answer  the  guestion; 
is  d^  (z,  P),<djj  ? If  this  can  be  answered  aff  irna  ti  vel y we  will 
find  this  out  at  the  point  where  any  upper  bound  is  less  than  . 
If  the  answer  is  no,  then  this  will  be  discovered  when  the  tree 
has  been  explored  to  the  point  where  all  itiniaun!  bounds  are 
greater  than  d . Similar  arguments  can  be  made  for  the 
qualitative  leve l-o f-e f f ort  required  to  answer:  is  d^(z,P)>dp? 
From  this  discussion  we  can  see  that  the  search  will  be 
inefficient  onlv  if  dpC^  dj(2,S(T))  and  a large  number  of  the 
strips  are  nearly  d^  from  z.  Figure  5a  shows  this  case  together 
wit-h  a more  representative  example. 


Figure  5.  Two  of  many  Possible  Geometries  When  Testing 
the  Distance  of  a Point  from  an  P. 

To  summarize  this  discission,  we  provide  the  algorithms  to  test 
for  d^(z,F)  < dp  and  d^{z,?)  > d^  . These  algorithms  use  the 
notion  of  the  distance  of  a point  to  a set  which  is  defined  as 
follows.  For  any  strip  S,  if  a point  is  outside  S i.e.  S 
then  its  distance  to  S is  characterized  bv  the  set  theoretic 
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distance  de(z,S)  = ain  d(x,z)  where  d is  the  euclidean  distacce 

xtS 

between  the  points  x and  z.  For  clarity,  the  algorithms  are 
presented  as  procedures  in  a pseudo-Algol  language.  Rigor  has 
been  sacrificed  mainly  in  the  specification  of  data  types,  tut 
these  should  be  obvious  from  the  earlier  definitions. 

Algorithm  A1:  Is  a point  within  dC  of  a polyline? 
boolean  procedure  Within  (z,dQ,T) 
begin 

if  dG<  ds(z,S('I))  + Z.deltaCT)  then  return  (true); 
if  z / S (T)  and  d1  >ds  (z  ,S  (T) ) then  return  (false); 
return  (Within  (z,  dO , L Son  (T)  ) or  Within  ( z ,d  C , F Son  (1)  ) ) ; 
end  ; 

Algorithm  A2:  Is  a point  further  than  dl  from  a polyline? 
boolean  procedure  Further  (z,dC,T) 
begin 

if  dC<  ds(z,S(T))  + 2. delta  (T)  then  return  (false); 
if  z ^ S(T)  and  dO>ds  ( z , S (T)  ) then  return  (true); 
return  (Further  (z,d  ) , LSon  (1)  ) and  Further  (z,  dO,  ESon  (1)  ) ) ; 
end ; 


3,2  Displaying  a Polyline  at  Different  Resolutions 

As  previously  demonstrated  in  Section  2,  a polyline  may  be 
represented  as  a set  of  strip  segments  such  that  each  strip 
segment-  L has  a resolution  delta  less  than  some  fixed  deltaO. 
The  algorithm  to  display  such  a representation  using  the  strip 
tree  is  as  follows.  This  algorithm  uses  a device-dependent 


r 
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5ut  rou*:  ine  Displavl-faC'*: angle  which  pain^?  the  rectangle  on  the 
par^ic'ilar  display  device, 

Alnorithni  .>3:  display  a polvline  at  Resolution  deltaO 
procedure  PolvDispiav  (TfdeltaO) 
begin 

if  delta  (T)  < deltaO  thtai  D ispla  yl  ec  ta  ng  le  ( L (T)  , delta  (I)  ) 

else  (PolyDisplay  (LSon  (T)  , deltaO)  and  PolyDisplay 
Son  (T)  , deltaO)  ) ; 
end  ; 


3.3  Intersecting  Two  Polylines 

C ne  of  the  important  features  of  the  representation  is  the 
ability  to  compute  intersections  between  polylines.  Strip  trees 
provide  the  facility  to  not  only  compute  intersection  points, 
tu*-,  in  he  case  where  lower  resolution  is  satisfactory,  to 
compute  small  areas  con'^aiaing  the  intersection  points  at  great 
computational  savings.  In  order  to  develop  the  intersection 
methodology,  we  need  the  following  definitions; 

A.  Twe  strip  segments  (LI  derived  from  ?1)  and 
(L2  derived  from  P2)  do  not  intersect  iff 
Ll/1  LL=  0 

E.  Twe  strip  segments  Li,  L2  have  a clear 
intersection  iff:11+  and  11-  intersect  12+  and 


12-. 
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C.  Two  strip  se^iuents  LI  and*  L2  have  a possible 
intersection  if  condition  D is  not  satisfied 
yet  Li  n 

These  cases  are  iilnstrated  tv  Fiqure  6.  fairly  obvious  tut 

very  important  lemma  is: 

Clear  Intersection  Lemma.  f-cucker,  1976]  If  two 
E-^rip  segments  have  a clear  intersection  and  the 
strips  are  both  compact,  then  the  corresponding  ?s 
must  also  intersect. 


To  see  this  for  condition  3,  consult  Figure  6b-  ?1  divides  the 
reqion  T into  two  parts  and  P2  must  cross  from  one  to  the  other. 
Th*-  onlv  wav  the  ?2  can  do  this  is  by  intersecting  PI. 


Figure  6:  Different  ways  Strips  can  Intersect 

The  algorithms  to  check  fer  intersections  between  two 
polvlines  are  recursive,  and  assume  the  existence  of  an  integer 
procedure  St  r ipl  n ^ ersec  * ion  which  will  return  the  type  of 
inte rs-'c tion  and,  in  the  case  of  a clear  type,  will  return  a 
parallelogram  Q containing  '■he  intersection  points. 


\lgori'-hm  \4: 


Finding  out 


w he  ‘■her 


wo  polylines 
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in’-crsfect 

Cr'.'ii tne nt . If  the  two  root  strip  segments  do  not 
u in*ersect  ♦■.hen  tht  Ps  do  not  intersect.  If  the 

roo*-  seqneiits  have  a clear  intersection  then  the 
?s  intersec*.  Since  the  task  is  to  just  determine 
wiisthor  or  not  an  intersection  exists,  we  are  done 
the  moment  we  find  a clear  intersection. 

I 

boolean  procedure  In'ersection  (II, 12,  Primitive 
Flaq) 

commen*  ?rimi*ive  Fiaq  allows  the  use  of  a single 
strip  as  the  first  argument 
begin 

Case  Striplntersection  (S  (T  1)  , S (T 2)  , Q)  into 
fNulll  return  (false), 

r Possible  1 if  ( \ r ea  (T 1 ) > J r ea  (12) ) or  (Primitive 
Fiaq)  then 

return  ( (Intersection  (L£on  (11)  ,”*2)  or 
(■ntersection(FSon(11)  ,12))  ; 
else  return 

(Intcrsection(Ii,LSori(T2))  or 
Inte  rsectior.  (1 1 ,H3on  (12)  ) ) ; 

[Clear!  re*urn  (true)  ; 
end ; 

This  orocedure  is  easily  modified  to  return  a set  of 
para lleloqrans  comprisinj  intersection  points. 


Further  easy 
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nod  i f ic.i t ionr  can  Le  .uade  *0  cousttain  these  paralleloqraaiE  tc  be 
of  a certain  size  related  to  the  delta  (Tl)  and  delta  (T2)  ; i-e., 

the  7 can  be  made  to  be  as  small  as  we  want. 

Kotr,  however,  tnat  smaller  resolutions  may  • be  much  icre 
coiunutat  icnally  expensive,  as  shown  in  the  following  example 
(Figure  “)  where  intersection  at  the  coarsest  resolution  is 
simple,  hut  niultiple  intcrcecticns  occur  at  lower  levels. 


Fi  rur  s : \r.  intersection  aav  he  simple  at  one  le''el  and  complicated 

at  lower  levels. 

If  two  Ps  are  not  convoluted  about  each  other  the 

intersection  will  bo  computed  in  0(mlog(n))  steps  where  m is  the 
numter  of  intersection  points.  If  the  ?s  do  not  intersect  but 
have  a closest  distance  d =ds(Pl,P2)  then  this  will  be  discovered 
it  1 l.vel  in  the  tree  no  deeper  than  a point  where  c)/m  4 deltal  + deltaZ. . 

The  worst  case  performance  is  intolerable  as  the  algorithm's 
v:ofii pu ti t ion  will  grow  exponentially  as  long  as  all  the  strip 
segments  in  one  t^^e  intersect  ail  the  strip  segments  in  the 
other.  In  fact,  .jg  computation  can  be  shown  to  be  0 (2  ) where 

" is  thr  sum  ot  the  depths  in  each  tree  where  the  comparisons  are 
t.ikinq  Dlice!  If  ‘his  situation  were  incoun+’orcd  in  a practical 
aorliciticn,  one  wav  of  handling  it  would  be  to  report  the 
njssitlT  inters^^c*  ion  rajions  at  the  point  where  ‘he  limit  of 
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some  Louiid  on  aiio^.tGi  resources  was  exceeded. 


3.4  Ihc-  Union  of  "^wo  Poiviines 

The-  union  of  rwo  strip  trees  can  be  accomplished  by  aefininq 
a s^rip  that  covers  both  of  the  two  root  strips. 

■ilqorithm  A5;  P-P  Uni  an. 

For  two  Ps  defined  tv  ^ . . . y^'J,  ^ 

^rcat  these  as  two  subsets  and  concatenate 
the  subsets.  That  is,  the  resultant  ordering 
is  such  that  we  have  v =y ' , v =v  ".  Mow 

o 0 *n4n4l  ^ 

define  a strip  seqment  that  covers  ‘ly^  7 • • • 


.7  ( such 

that  c= 

0 and  delta 

• 

II 

ny 

corjS*ruc  tion , 

this 

satisfies 

a 1 i 

the 

properties  of  a strip  sequent.  Maxe  this  the 
root  node  of  a new  p-tree.  The  ♦wo  subtrees 
are  ♦he  ♦wo  ?s  of  the  union. 

This  construction  is  shown  in  Fiqure  The  variable  c is 
defined  below. 


Fiqure  8:  Const ruct  ion  for  Unior.  of  Strip  Traes 


Pepresentinq  Two  Polylines 
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Cf  cour5=€  ♦■his  constcuct ion  introduces  a problaa  in  that  the  new 
strip  is  no  lonqer  compact  and  therefore  the  Clear  Intersection 
lemma  no  lonqer  holds.  To  overcome  this  problem  we  must  add  one 
bit  of  information  to  each  node  to  nark  whether  the  underlyinq 
polyline  is  compact.  Since  later  algorithms  may  result  in 
iinderlvinq  polylines  that  are  disconnected,  we  include  this  in 
♦•he  follow ir.q  definition  of  C: 

r ('")  =(  1 P represented  by  S (T)  is  known  to  be  compact  and 

J connected 
I 0 otherwise 

W i*  h this  strateqy  we  can  preserve  the  eloquence  of  the  previous 
alqori-i-hns  in  ♦•lie  folj.owinq  manner.  lnhen  bit  C (T)  is  not  one  we 
apniy  ♦•he  recursion  rejardless  of  the  intersection  type.  In 
alqorirhm  iU  tliis  means  that  clear  intersections  are  reported  as 
possil  le  if  ♦•he  bit  C(T)  is  set. 

This  technique  ■:aa  also  he  used  as  a diqitization  method  for 
□ non-ccnnected  segments  , 'd,');  Cy;,,  , • , 

These  scqments  are  qiven  an  orderinq  as  shown.  The  previous 
d iait  i?  a*  ion  algorithm  is  applied  to  this  set  af  points,  and  the 
perpendicular  dis<-ance  d^  is  computed  from  the  set  of 
disconnected  vs  ani  used  to  define  the  of  the  root  strip  as 
befare.  Pj-wcver  now  t.ie  set  is  divided  into  two  subsets  of 
connect  .d  segments  (rather  ♦•han  using  v*)  and  the  d iqit iza ♦: ion 
ai?ori‘hm  is  applied  recursively  to  the  subsets.  Cnee  this 
proce.'s  produces  connecfad  subsets,  the  earlier  diqitization 
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schcnt  i£  applied. 

4.  Areas  FepreseLted  by  Strip  Trees 

We  ♦•ake  the  boundary  oi  an  area  to  be  a closed  polyline, 
-at eres* inqly  encuqh,  the  digitization  method  described  in 
Section  2 works  for  closed  polylines  and,  incidentally,  also  for 
sfclt-ir.tersectinq  polylines.  Furthermore,  if  an  area  is  not 
simply  connected  it  can  still  be  represented  as  a strip  tree, 
which  a*-  some  level  has  connected  primitives.  The  method  for 
doing  so  was  described  in  the  previous  section.  If  a reqion  has 
holes  it  can  be  represented  by  a single  boundary  curve  using  a 
construction  ^Figure 


Figure  1;  \ region  with  a Hole 

If  *hc  holes  are  inportant,  they  themselves  should  be 
independently  represen*ei  as  strip  trees. 

'^he  most  remark  sole  fact  is  tha*-  by  representing  an  area  in 
*-his  way  many  useful  operations  such  as  intersection  between  a 
polvline  and  an  area,  ietirmining  wnether  a point  is  inside  an 
trea , and  intcrsectiui  two  ircas  are  carried  cut  vsrv 
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*:  f f icic  i;-^  iy  . 

U.1  Det'irm  ininq  Whether  i Point  is  Inside  an  Area 

The  strip  tree  representation  of  an  area  by  its  boundary 
allows  * he  deter  a inat ion  of  whether  a point  is  inside  the  area  in 
a straiqhtf erward  oiaaner.  If  any  semi-infinite  line  terminatinq 
at  the  point  intersects  the  boundary  of  the  area  an  odd  number  of 
times,  the  point  is  inside.  This  result  appears  in  f Minsky  and 
Papert,  19f:91.  This  result  is  computationally  simplified  for 
strip  trees  in  the  followinq  manner; 

Point  Membershio  Property 

To  decide  whether  a point  z is  member  of  an 
area  repres^intod  by  a strip  tree,  we  need 
only  compute  tae  number  of  c^ear 
intersect  ion  i of  the  strip  tree  with  any 
semi-infinite  strip  I which  has  delta  = 0 and 
emana+:os  from  z.  If  this  number  is  odd  then 
the  point  is  inside  the  area. 

An  extension  to  the  clear  intersection  lemma  which  makes 
this  property  hold  is  that  the  underlvinq  curves  may  intersect 
more  than  once  tut  must  intersect  an  odd  number  of  times.  The 
rollowinq  aljorithm  is  used  to  determine  whether  a point  is 
inside  an  a re  a; 

Alqori‘l;in  A5:  ?oin^  Member  shin 


▲ 
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booL'^ati  proctdur‘1  Inside  (z,T) 

1 7 i I. 

Crc-ateS*:rip(3Q,z) 

CDiniocRt  'Zri.a  t cSt  r ip  creates  a strip  for  the  half  line, 
id  N cC  f V learl  nter.SteCt  ions  (5  C , 1)  is  odd  then  return  (true) 
else  return  (fais-=>)  ; 

end;  f 

intejer  procedure  NoOf  Cie arin t ersect  ions  (S, '^) 
fceqin 

Ca  se  S*-. ripint  G rsection  (S  ,S  CJ)  ) in^  o 

1 

[ i:u  11  ] ret  urn  (0  ) ; J 

rPossihle)  return  ( IloC  fClearl  n ter  sect  ions  (S  ,LS'on  (T ) ) , 

i 

♦ NoCf Clearlnterscctions  ( S, ” Son  (1) ) ) ; 
r Clear  1 return  ( 1 ) ; 
e nd  ; 

A potential  difficulty  exists  with  the  procedure 
;ior  f Clear Interse  G+ ions  when  tj^t  strip  SO  is  tanqcnt  to  the 
polviine.  Since  ^his  oroblem  will  only  occur  at  the  lowest  level 
of  ‘he  tree,  .e  can  exanine  neiqhbcrinij  leaves  of  the  tree  to 
rerolvt  it. 


h.2  Intr.  rsectinq  a polviine  with  an  Area 


Che 

poi V 1 int 
new  Tee 
'-his  can 


St  rate  qv 
w i ♦:  h a 
f f ) r the 
I e 


behind  intersfcctinq  a strip  tree  representir.q  a 
strin  tree  representing  an  area  is  to  create  a 
porti  in  of  the  polyline  which  overlaps  the  area, 
‘rimminq  ‘be  oriqinai  polyline  strip  tree. 


done-  L V 
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“tiis  i'  Oone  efficiently  Ly  ♦nkir.q  advantage  of  in  obvious 
proper'-v  of  ♦ho  interdiction  process: 

Pruning  ?rop:-rty:  Consider  two  s-^.rips  Sp  e 


Tp 

and  Sa  e Ta. 

J.f  the 

3p  0 

Ta  is  null,  then 

(a) 

if  any  point 

on  Sp 

is 

inside  Ta  the 

en*- 

ire  t ree  w hos 

c roo  t 

strip 

is  Sp  IS  inside 

or  on  Ta  and  (b)  if  any  point  on  £p  is 
outside  of  Ti  then  the  entire  tree  whcse  root 
s*-rip  is  £p  is  outside  of  Ta . 

This  leads  to  the  recursive  procedure  \7  for  polyline-area 
intersection  using  trees.  Mote  that  since  strip  nodes  under  a 
clear  or  possible  strip  intersection  may  be  pruned,  the  bi*  c for 
•^he  la- ter  strip  is  sot  to  0 to  denote  that  it  no  longer  has  the 
coQoac*ness  property.  Of  course  as  repeated  intersections  are 
carried  out  with  different  areas  more  and  more  upper- level  strips 
mav  have  their  bits  set  to  0;  nevertheless,  the  intersected 
polyline  is  accurately  represen'^ed  at  ‘•he  leaves  of  the  strip 
“ r ee . 

Note  ‘■hat  if  the  ooivline  strip  is  "fa‘-ter,"  i.e.  , Area(TI) 
> Area(T2)  , wc  car.  copy  the  node  and  resolve  the  intersection  at 
lower  levels,  whereas  in  the  converse  case  we  have  to 
.-eguentia  lly  prune  tne  tree  by  first  intersecting  the  polyline 
strip  with  the  left  area  s^rip  and  then  intersecting  the 
rc.su  Ltant  pruned  tree  witii  the  right  area  strip. 


.iiiIori*l..Ti  A":  Foi  vl  inc- lut  rr  st-c*- ior; 
r et *-*reiice;  proci-aar^  ? civ 'i  r-sal  n*  (1  1 , 1 2) 
b a i n 
A : =1  2 

coimp.cn*  A is  a cjiobai  used  bv  FAIn^; 
return  (PAIr.t  ("'1, 12)  : 
eni  : 

retereiicfc  procedure  ? ilnr  (T  1 ,1  2) 
be  qiri 

Case  Stripiri*  (11  ,T2)  into 
fliull  nr  Frimitivi] 

it  latersectioa  (11, A,  IP CS)  = null 
if  Insi  ds  (1  1 , A)  then  return  (1 
eire  return  (null)  ; 
e Lsc  return  (11); 

fZl'^ar  or  Forsibial  if  Area  (1  1)  > Area  (12) 
t € 0 i n 
C (f:i)  : = 0 

cenreent  iion-coapact  strip 


XEeq  (III)  : = 

X2e  ({ID; 

YPeq  (11)  : = 

Y3o  ) (1 1)  ; 

X Zn  i C-I'")  : = 

XZnl  (11)  ; 

YZnd  (HI)  : = 

Yind  C"  1)  ; 

Are.1  (i:i)  : = 

Area  (1  1)  ; 

ISon  (’;i)  ; = 

P xint  (LSon  ( :i)  ,12)  : 

' Son  (111  : = 

Flint  (-3on  (1  1)  ,T2)  : 

t hen 
1) 

t hen 
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t Ice  GCBJicPi*  ir9a(T1)<  Area  (T  2) 

! fc^urn  (’’AInt  (?\Int  (T  1,  LSoii  (T2)  ) , ?S^n  (12)  ) ) ; 


end; 


4.  A In r erstct inq  Iwo  Areas 

1!i9  prctlcm  of  in ‘■.ersectinq  f «c  areas  can  he  efficientiv 
carried  oj*  usinq  thair  srrip  tree  representations.  The  method 
is  to  decompose  the  prohlem  into  two  polyline  area  intersection 
prohlems  (refer  to  Piqare  10). 


’^’ij'ire  10;  Decomposition  of  Arei-Area  Intersections 

If  wc  treat  the  hoaniiarv  5f  A1  as  reorssentinq  a polyline  instead 
of  re  or  cs’^r.*- inq  an  area  and  intersect  its  strip  tree  with  the 
strip  ‘ree  re- prese  ntin  j A2  tne  lowest  level  result  is  shewn  hv 


* he 

h i c k 

line 

in  Fi  ju-'- 

1 Oa  . 

If  we 

reverse  the  roles 

of  the 

t wo 

St 

rip 

*■  reen 

*he  result  is 

qiven  hy  t 

he  thick  lines  in 

Fi q ure 

1 01-.  - 

The 

u nion 

o^ 

these  -.wo 

strip  trees 

(soe  Sect  ion  3.  4) 

is  the 

a ns  w 

'.r 

w fc 

want! 

h u s w £ ca  n 

write 

t he 

area-area  inter 

section 

proc 

eu 

urc 

in  terms 

of  strips 

as  follows; 

Aljori'hm  A’’:  A.rra-Arid  Intersection 


rrfercr.cc  procedure  AreiArealn*:  (Tl,""!) 


A 
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I eg  in 

return  (Ir.ion  (PoivArealar  (T  1,?2)  ) , (Poly  .\realnt 
r2,Z^)))^, 
e-nd  : 


whtre  Jiiicn  is  a procedure  that  acco mplislie s the  construction 
described  in  Sec'-ion  3.4. 


Note  that  in  ■'■he  case  of  areas  that 
fraa-Tifeiits  their  Loundaries,  the  erder  of 
preserved  by  the  intersection  procedure, 
were  auarantecd  '•hat  strips  in  the. 
accordinq  tc  the  arc  lonqth  of  their 
'dowever,  all  the  ether  properties  of 


intersect  in  a wav  that 
the  segments  will  not  he 
('In  til  this  pcint  we 
tree  would  be  ordered 
underlying  polylines), 
the  re- presen  tation  are 


pre  se  r ved  . 


4.4  "ho  Union  Operatijn 

The  union  operations  are  sliqh‘ly  simpler  than  the 
intersection  operated'..  For  the  union  of  a E-tree  and  a P-'Tee 
we  use  a construction  similar  to  "ho  digitization  methods  for 
disconnc*cd  Ps.  The  rosult  is  a ?-*refe.  Note  that  the  union 
operation  fci  strip  treos  is  not  commutative.  \lso,  we  do  not 
lefine  a union  ofierati-jn  for  a strip  "tee  representing  a polyline 
and  a -"rip  "ree  rcpresen*ing  a region.  The  union  of  two  region 
strip  tices  is  definei  and  is  a region  3"rip  tree.  Tf  these  two 
strip  "tees  do  no"  int-3rsect,  then  the  union  is  s"raiqht  ror  war  d 
am;  is  idcn"ical  to  t h«  met  nod  for  polylines.  Uowever,  if  the 
cor.trarv  is  true,  then  we  miist  go  to  "he  t^^ubie  ot  itfining  a 
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new  strip  tree  -hat  reprasents  the  anion  by  finding  the  points  of 
intersection  in  the  saae  way  as  was  done  for  region  strip  tre-e 
in- e rsect ions. 


I 
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' 5.  ConcluFions 

I 

I Strip  trets  provide  i powerful  represent ar  ion  for  polylines 

y 

^ and  areas.  Current  work  is  directed  towards  characterizing  their 

conputational  complexity  more  precisely  but  it  can  already  be 
shown  that  the  representation  is  superior  to  its  competitors. 
The  main  drawback  is  that  there  is  a large  overhead  in  terms  of 
space.  If  n is  the  required  space  to  represent  a polyline  then 
its  strip  tree  will  taxs  about  4n  space  units.  Mso  the  creation 
of  a s-^rip  tree  is  a laborious  process,  requiring  0 {n  log  l)  tine 
units.  However,  neither  of  -^hese  drawbacks  are  thought  to  be 
important  in  the  use  of  this  representation  for  geographical  data 
bases . 

The  representation  defines  strip  segments  as  primitives  to 
cover  subsets  of  the  line  after  fTeucker,  IP'ol.  Cur 
or ganira+ ion  of  -^hese  segments  into  a tree  may  be  viewed  as  a 
particular  case  of  a general  strategy  of  dividing  features  up  and 
covering  them  with  arbitrary  shapes  such  as  depicted  by  Figure 
10.  ether  attempts  in  this  class  have  been  tried  by  f Barrow  et 
ai.,1'^T"';  Burton,  19^7;  Tanimoto,  1 975  1,  but  thev  do  not 
cao^-ure  the  notions  of  orientation  and  resolution  anywhere  nearly 
as  precisely  as  strip  segments,  and  do  not  have  the  union  and 
inte rsec icn  properties. 
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Figure  4.  A portion  of  an  encod 


Figure  5: 
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Figure 
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Two  of  Many  Possible  Geometries  When  Testing 
the  Distance  of  a Point  from  a P. 


6:  Different  Ways  Strips  Can  Intersect 


I 


nple  at  one  level 
levels. 


on  of  Strip  Trees 
:nes 


a Hole 


